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Abstract 

In highly correlated systems one can define an optical self energy in analogy to its quasiparticle 
(QP) self energy counterpart. This quantity provides useful information on the nature of the 
excitations involved in inelastic scattering processes. Here we calculate the self energy of the 
intraband optical transitions in graphene originating in the electron-electron interaction (EEI) as 
well as electron-phonon interaction (EPI). Although optics involves an average over all momenta 
[k) of the charge carriers, the structure in the optical self energy is nevertheless found to mirror 
mainly that of the corresponding quasiparticles for k equal to or near the Fermi momentum kp- 
Consequently plasmaronic structures which are associated with momenta near the Dirac point at 

= are not important in the intraband optical response. While the structure of the electron- 
phonon interaction (EPI) refiects the sharp peaks of the phonon density of states, the excitation 
spectrum associated with the electron-electron interaction is in comparison structureless and fiat 
and extends over an energy range which scales linearly with the value of the chemical potential. 
Modulations seen on the edge of the interband optical conductivity as it rises towards its universal 
background value are traced to structure in the quasiparticle self energies around kp oi the lower 
Dirac cone associated with the occupied states. 
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I. INTRODUCTION 



The isolation of a single layer of graphend^ lead directly to the establishment of the 
charge carriers in this system as effectively massless Dirac fermions. The carrier mobility in 
graphene is particularly large and the number of carriers can be modified through the gate 
voltage in a field effect configuration. It is also considered promising for the development of 
a new generation of electronic device^. Several reviewd^lS! have appeared which document 
the unusual charge transport in this material. 

The optical properties of graphene are particularly interesting and an experimental review 
was given by Orlita and PotemskP. To summarize, graphene shows a Drude-like intraband 
conductivity centered around i7 = which is followed at higher energies by a region of 
reduced conductivity before the sharp rise to a constant universal background value of 
o"o = vre^/ {2hj^^^^ (e is the charge on the electron and h Planck's constant) at Q equal 
to twice the value of the chemical potential /i. This part of graphene's optical response is 
due to the interband transitions between the fully occupied lower Dirac cone of the valence 
band and the unoccupied part of the conduction band Dirac cone. Another interesting 
optical property is a giant Faraday rotation. As demonstrated by Crassee et alJ ^^ * ^^ *, a single 
layer of graphene turns the polarization of incident light by several degrees in relatively 
modest magnetic fields. This property could be useful for ultra thin infrared magneto-optical 
devices. A related effect is the enhanced dichroism of graphene nanoribbon^ where the 
polarization state of the light effects its absorption. While the effects described above, and 
many others, can be understood within a bare band picture, some cannot. As an example, 
a large amount of absorption, of the order of one third the universal background value 
(To = vre^/ {2h) is observed in the Pauli blocked region of the spectrum which lies between 
the intraband Drude peak and the universal background which occurs above Q > 2fi due to 
interband transitions. This additional absorption can be understood to arise partially from 
the electron-phonon interactions (EFVp^ with an additional contribution from electron- 
electron (EEI) correlation^. These many-body effects provide a finite self energy part to 
the electronic spectral density, A{k,u), for momentum, k, and energy ,u. This quantity 
is measured directly in angular resolved photoemission experiments (ARPES)P2H25! the 
occupied values of k. Recent ARPES data on graphene has revealed the splitting of the 
Dirac point into two separate points with a region in between showing plasmaronic features 
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associated with a coupling of electronic states with plasmons. Such structures are also 
predicted to exist in the density of states, N^u), which depends on an average of A{k,u) 
over all momentum values, k, but have yet to be clearly isolated from other many body 
effects in scanning tunnelling microscopy (STM)?^'^. In this regard graphene is different 
from a conventional metal in which the renormalized density of states is not expected to show 
sharp many body signatures, in as much as its bare band value does not vary significantly on 
the energy scale of the exchanged boson involved-l^^^^^owever, the linear in 00 dependence 
of N{u) in graphene is sufficient to modify this expectation. 

Of course there are many known effects of electron-electron (EEI) and electron-phonon 
(EPI) interactions which can provide profound modifications of the physical properties of 
metals. For example they lead to inelastic scatterin^^SES g^^^j ^]-^g superconducting state 
they provide the so called strong coupling correctiond^lHlIl canonical Bardeen-Cooper- 
Schrieffer (BCS) theory. These corrections are remarkably well described by Eliashberg 
theory for electron-phonon coupling. Inelastic scattering is quite distinct from other possible 
complications such as the presence of van-Hove singularitie^^^BS ^nd even strong scattering 
anisotropieJ^^'^ since these complexities are not expected to be important in discussion of 
the low energy properties of graphene. 

In this paper we calculated the effect of electron-electron interactions on the optical self 
energ}^ of the massless Dirac fermions of graphene. The optical self energy is defined 
in direct analogy to the quasiparticle self energy. It is however a two body property and 
is defined in terms of a complex generalized Drude form in which both carrier effective 
mass and scattering rate acquire a temperature and frequency dependence. Optical mass 
renormalization and scattering rate are the central elements to be calculated here for the 
intraband transitions. For comparison with the EEI case, and to provide additional insight 
into that case, we also study the effect of the EPI on the same quantities. This involves 
the exchange of phonons and provides well known many-body corrections to Drude theory 
in conventional metals, such as the Holstein phonon assisted absorption sidebandd^^'^. In 
this case details of the electron-phonon spectral density denoted by a^F^u) are encoded in 
the optical mass and scattering rate and these have been "inverted"'^ to yield a'^F{uj) which 
is a dimensionless function closely related to the phonon distribution F{uj). The difference 
between these two fundamental functions lies in that each phonon in a'^F{uj) is further 
weighted by an appropriate factor related to the electron-phonon interaction strengthP^. 
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By contrast, for EEI it is not from the outset guaranteed that such a boson exchange concept 
can be apphed to the excitations involved, be they particle-hole, plasmons or plasmarons. 
Nevertheless it will turn out to be helpful in the interpretation of the structure seen in 
optical properties when the electron-electron interactions are included, to understand how 
such a concept can be approximately applied to the results. 

Along with the intraband optical transitions which provide a Drude-like respond center 
around uj — and which are analyzed in detail here, graphene also exhibits interband 
transitions. These set in only at photon energies, uu, comparable to twice the value of the 
chemical potential, 2/x, and provide a nearly constant universal background uo = Tre^/{2h). 
These optical transitions are not part of the usual single-band metallic case and cannot be 
naturally analyzed with an extended Drude form. Nevertheless we find that the structure 
seen in the low energy rising edge of the interband piece can also be related to well defined 
excitations of the electrons. 

In section II we define the optical self energy, T,°p{uj), and relate it to the optical con- 
ductivity which can be calculated from its Kubo formula as an appropriate overlap of two 
electronic spectral functions for the same momentum but with energy arguments displaced by 
the photon energy, ou. We give the formulas needed to calculate the quasiparticle self energy, 
'Ef{k,uj), which enters the spectral function. In general T:f{k,uj) depends on momentum, 
k, energy, ou, and on valley index, s — +/— for conduction and valence band respectively. 
To include electron-electron interactions we employ a random phase approximation for the 
dynamical screening of the potential between electrons. This latter quantity is reduced by 
the medium dielectric constant which enters in its denominator as an average value of the 
materials above and below the graphene sheet. We present our numerical results for pure 
EEI, EPI and combined EEI-I-EPI, and begin by emphasizing the structure that enters the 
optical mass renormalization, A°^(a;), as a function of u. This structure reflects details of 
the excitation spectrum to which the massless Dirac electrons are coupled. We compare 
the structure seen in X°p{u) with the structure that enters its counterpart, X'^^\k,uj), which 
can depend on momentum as well as on u and emphasize that only its value at or near 
k — kp appears to enter importantly into y^iuj). In section III we introduce simplified, 
approximate but analytic formulas which apply to the case of a general energy dependent 
electronic density of state and specialize them to the specific case of graphene. We obtain 
formulas for real and imaginary part of the quasiparticle and optical self energy in a gen- 
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eral boson exchange model and compare them. Analytic results are given for their zero 
frequency limits {u = 0). Section IV deals with the problems of fitting quasiparticle and 
optical self energy data obtained from complete calculations based on the Kubo formula, to 
our simplified equations for a boson exchange mechanism. A maximum entropy inversion 
technique is used to obtain the best fitting electron-boson spectral density. The consistency 
of such an approach is checked, and its limitations commented on. Despite these limitation 
the procedure yields insight into the nature of the excitation spectrum that provides the 
scattering when electron-electron correlations are included. In section V we consider the 
case of the boson structure that arises in the interband part of the optical conductivity of 
graphene. Section VI contains our conclusions. 



II. INTRABAND OPTICAL SELF ENERGY 

It has become common to analyze the optical conductivity, a(T,u), of a correlated elec- 
tron system in terms of a generalized Drude form in which an optical self energy, S°^(T, u), 
is introduced. Here, T is temperature and u is photon energy. We writ^ 

a{T,u) ^ a,{T,u) + ^a,{T,u) = ^^J^^^^^ (1) 

where ai{T,u) and a2{T,u) are the real and imaginary parts of the conductivity and Qp is 
the plasma energy. It is further instructive to write E°^(T, u) = E^^ + iEg^ in terms of its real 
and imaginary part with — 2S2^(T, w) = l/r°^(T, w) and — 2S°^(T, u;) = u[m*p{T,uj)/'m ~ 1] 
where t°'p{T,uj) is an optical scattering time and m*^p{T , uj) j m — 1 = A°^(T, w) is an optical 
mass renormalization factor. In terms of these quantities the conductivity takes on its usual 
non-interacting form but now with frequency dependent scattering rate and effective mass. 
In particular, the real part of a{T, uj) is 

,(rn\ = ^ (0\ 

47r[a;(l + A°P(T,a;))]2 + (l/r°P(T,u;))2- ^> 

For graphene the Kubo formula gives the conductivity of the massless Dirac fermions in 
terms of their spectral density, A^{k,uj)^ where k is momentum and u is energy. The 
'±' notation refers to the conduction and valence bands respectively. For the bare bands, 



e±{k) = ±VF\k\ with vp the Fermi velocity. These bands define two Dirac cones which meet 
at k = 0, the Dirac point. In graphene there are two such points in the Brillouin zone at K 
and K' points which provide a valley degeneracy factor, = 2, in addition to the standard 
spin degeneracy (?s = 2.'^ 

The the real part of the conductivity is given by 



ai{T,oj) 4 



/oo rWc 
duj'[f{uj')~ f{uj' + uj)] / kdkA{k,uj')A{k,uj' + uj), (3) 
-oo Jo 



where Wc is the energy of the band cutoff. The imaginary part of the conductivity can be 
obtained by using a Kramers-Kronig relation. 

(T2{T,u) = P — -du (4) 

where P stands for the principal value integration. 

Here, the total charge carrier spectral function, A{k,uj) is given bj^ 

Mku) = S-^- -i^^Tik^^) 

ti'^i^- ReSf (A;, u) - e^]^ + [ImSf (A;, u)f ^ > 

where = svF\k\ — /io with vp the Fermi velocity, /io the chemical potential and T,f{k,uj) 
the quasiparticle self energy due to many body interactions which can depend on the band 
index, s. Here, for simplicity, we will be interested in treating the case of an electron- 
phonon interaction with coupling to a single Einstein mode at energy u = 00 e^^^ Detailed 
calculations of the electronic self energies in graphene due to electron-phonon coupling have 
been done in density functional theory.'^^^ An important observation made in referencd^ 
for the present work is that the results of such complex computations have little dependence 
on the direction and magnitude of the electron momentum k and can be modeled in a 
first approximation through coupling to a single phonon mode at energy u = 200 meV. 
Here we adopt this model but will also allow for coupling (still assumed independent of 
direction and magnitude of k) to a group of phonons rather than a single mode. This can 
be accomplished by the introduction of an electron-boson spectral density a'^F{u) again 
assumed to be independent of the Dirac fermion momentum k. With such a simplified 
model, the electron self energy is momentum and band index independent and follows from 



the equation!^ 

Jo J^oo 

with N{ll!')/No the density of electronic states which for graphene we take as \u' + /io|//^o 
as a first approximation and the Heaviside function <d{uj). There is no dependence of this 
self energy on the electron momentum, k. 



+ 



u' + iO+ uj + u -oj' + iiO+ 



(6) 



The electron-electron interaction provides a second contribution to the self energj^^^^ 
of the massless Dirac quasiparticles in graphene which is more complicated as it depends 
on electron momentum and band index in an essential way. However, there is a simplified 
scahng which applies. One can show that for reduced variables, k = k/kp, u = u/fio, 
Ils{k,u) is a unique function where = T.^/ fiQ. The expression for the self energy has two 
pieces and is written as 



s'=±l 



2n 



dqdO 



kq 



a 



27r g 



dn 

2^ 



k+q 



in 



^ k+q ' 



+ [e 



k+q 



-CO 



where 



-Fes' ( 



1 



'kk' 



[1 + cos{9^p)ss' 



(7) 



(8) 



and £ ^ is the inverse dielectric function. It is calculated here in a random phase approxi- 
mation and takes the form 



(9) 



q — all{q, n) 

with n(g, 0) the polarization function written in reduced variables.!^ The parameter a = 
ge^/leoVp) gives the overall strength of the coulomb interaction and is inversely proportional 
to the dielectric function of the substrate eo (average of top and bottom medium) and the 
degeneracy factor g = 4 = gsgv 



From Eq. ([T]) and the definitions of and we arrive at expressions for these two 
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optical constants which form the basis of our discussion, namely 



' %4^\ (10) 



We will apply these equations only for the intraband part of the conductivity of Eq. (|3]). 
From Eq. (jsj) we see that A(k, u) = A'^{k, u) + A~(k, u) where A'^ describes the upper Dirac 
cone and A~ the lower cone. The product of A{k, uj)A{k, u+Q) then expands into four terms, 
A'^ {k , u) A'^ {k , oj + Q) and A~{k,uj)A~{k,u + Q) the intraband terms and A'^{k,u)A~{k,u + 
Q) and A~{k,uj)A'^{k,u) + Q) are interband terms describing transitions from valence to 
conduction band. These will be discussed separately later. For now we deal only with the 
intraband part for which it is natural and useful to define an optical scattering rate and 
optical effective mass. 

In Fig. [l]we show results for — 2E°^(u;) (top row) and —2Tj2^(u) (bottom row) both in 
meV as a function of u over an small (left column) and expanded (right column) energy 
range (1 eV and 22 eV respectively) which extends well beyond the cut off of the bare band. 
Four distinct cases are shown. Pure EPI case, solid black line, the pure EEI case, dashed 
red line, and the results when the electronic self energy is fixed in momentum aX k = kp 
which we denote by EEIk=kp and is shown as the green dotted-dashed curve and finally, 
the EEI+EPI case shown as the double-dotted-dashed blue curve. In all cases the l/r°^(a;) 
remains positive definite for all u and is interpreted physically to be a scattering rate. By 
contrast, the real part of the optical self energy (— 2S°^(a;)) changes signs. In this paper we 
are mainly interested in the low energy part of this data, emphasized in the left column, 
rather than on how interactions modify the high energy band edge. It is important to 
keep in mind that for bare bands, the optical self energy, Zl°P(a;), would be zero. Therefore 
the deviations of the curves in Fig. [T] from zero carry the information on the effects of 
correlations. An interesting case to emphasize should be the elastic impurity scattering rate 
which gives the usual Drude theory. In this case the Kramers-Kronig transformation of a 
constant scattering gives a zero real part. 

In Fig. [2] we show results for the optical effective mass enhancement parameter, A°^((X'), 
derived from the data of Fig. [T] in the range up to 1.0 eV. We begin our discussion with 




FIG. 1. (Color online) The negative of twice the real part —2'E°^{uj) (top row) and imaginary 
part —2^2^(0;) (bottom row) of the optical self energy calculated from equations (10) and (11) for 
cr{Lo) given by the the Kubo formula, Eq. ([3]), for the real part of the conductivity and a Kramers- 
Kronig transform for the imaginary part. Various cases are considered. A pure electron-phonon 
case (solid black), a pure electron-electron case (dashed red), a combined EEI+EPI case (dashed 
double dotted blue) and a simulation of EEI with quasiparticle self energy approximated for all 
/c's by its value pinned at k = kp (dashed dotted green curve). The left column ranges up to 1 eV 
while the right column extends up to 22 eV. 
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FIG. 2. (Color online) The optical effective mass renormalization X°p{uj) vs oj for pure EPI (solid 
black), EEI alone (dashed red), combined EPI+EEI (dashed-double dotted blue) and a simulation 
of EEI alone for which the quasiparticle self energy is approximated for all k^s by its value pinned 
at fc = kp. 
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the EPI only case (solid black). To compute this curve we have used in Eq. ([6]) a delta 
function at w = fi^ to represent coupling to a single boson. Note the strong peak feature 
in X°P{u) at u = Qe- This quantity mirrors the sharp structure in the underlying spectral 
density, a'^F{u), and can be used as a first view of the boson exchange spectra involved in 
the quasiparticle scattering. A rule of thumb is that a peak in the electron-phonon spectral 
density translates into a corresponding peak in the effective mass renormalization X°p{u)) at 
this same energy. Also note the rather fiat value of below u = Qe with X"p{u) = 0.32 to 
be compared with a value of ~ 0.29 at energies above the phonon. This is to be contrasted 
with the case for the pure EEI, dashed-red, as well as the results of a simulation, dashed- 
dotted green, for which we ignore the variation of the coulomb self energy with momentum 
and have instead fixed it to its value at the Fermi momentum, k = kF- In comparison to 
the black curve, shows no peak but rather decreases smoothly from its value near a; = 0. 
This indicates that the effective boson spectrum to which the electrons are coupled is rather 
unstructured and fiat in its energy dependence. Including an additional phonon part to 
the electron-boson spectral density restores structure to the optical mass renormalization at 
UJ = Qe as we see in the dash-double-dotted blue curve. 

While we can make conclusions about the nature of the excitation spectrum involved in 
the optical renormalization from the resulting shape displayed by X°^{u) vs u we can get 
additional insight from the study of simulations based on a boson exchange model similar 
to the electron-phonon model but now for arbitrary choice of a'^F{u) chosen to simulate, as 
well as such a model can, the effects of the EEI. This will be comparable to previous results 
for marginal Fermi liquid systems. One can also gain insight by comparing our results 
for the mass renormalization A°^(a;) vs u obtained from optics with similar results for the 
quasiparticle. The quasiparticle self energy, T,'^P{k,u), depends both on momentum, k, as 
well as on energy, u, in contrast to its optical analog which depends only on energy since the 
conductivity of Eq. ^ involves an integration over k. Results for the quasiparticle effective 
mass X'^^{k, u) = ^[S'^*'(/c, u) — T,'^P{k, u = 0)] are shown in Fig. [sjfor four values of k namely 
k = kF (sohd black line), k = O.Gkp (dashed green line), k = O.Skp (dashed-dotted blue) 
and k = 0.0 (doubled dotted dashed red line). Note that the overall behavior of X'^^ik.u) 
for k = kF is very similar to that found in the pervious Fig. |2] dashed dotted green curve 
for the optical effective mass renormalization A°^(co'), which applies for a value of chemical 
potential equal to 400 meV. While the magnitude of the corresponding self energy is larger by 
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FIG. 3. (Color online) The quasiparticle electronic mass renormalization X''P{k,uj) for four values 
of momentum as a function of reduced energy u/ fiQ. The cases are k = kp (solid black), k = Q.Qkp 
(dashed green), k = O.S/cir (dash-dotted blue) and A; = (dashed double-dotted red). 



about 1/3, its variation with uj is much the same. More importantly as we move away from 
k = kp towards zero momentum, the corresponding structure in the quasiparticle curves 
does not look like the optical quantity. This is a first indication that optics (intraband) is 
most sensitive to the boson structure or excitation spectrum which dresses the quasiparticle 
at the Fermi surface, k = kp, rather than at the Dirac point near k = 0. This idea is 
investigated in the next section through maximum entropy inversion techniques. 



III. APPROXIMATE BOSON EXCHANGE MODEL 



To help us understand better the many-body interactions at play in graphene we have 
found it useful to introduce an approximate but analytic boson exchange model which was 
previously derived and applies for any energy dependent density of state N{uj). Here we 
specialize it to the specific case of graphene for which N{ijj)/No = \u + fio\/ fiQ. The formulas 
for quasiparticle self energy given by Eq. (|6| can be worked out to read 
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For the quasiparticle mass renormalization parameter, A'^^(aj), given by 



- wA^P(w) = ReS^P(w) - ReE'^f (0) (13) 
we evaluate the case of a; = to obtain 

X'^P^co = 0) = 2 rdna'F{n){]- - - + , + -In ]. (14) 

We emphasize that the first term gives the spectral effective mass A = 2 a^FiVt) /VLdVL 
which is the only quantity that would enter in a conventional metal with constant density 
of states. In general the quasiparticle mass renormalization, \^p{uj = 0), differs from its 
spectral value, X°^{uj = 0). The imaginary part of the quasiparticle self energy of Eq. g 
which gives the quasiparticle scattering rate can also be simplified. It reads 

- ImS'?P(w) = 71 dna^F{n)N{uj - n) (15) 
Jo 

for u > and 

-imj:ip{u) = TT [ dna^F{n)N{-\u\ + n) (le) 

Jo 

for a; < 0. Note that here we use the density of states of doped graphene with a linear band 
approximation as N^u) = |w + /io|//io- The resulting quasiparticle scattering rate is not 
necessarily the same for negative and positive u. In the limit of a; —t- we assume that the 
boson spectral density contains a constant u;° term, a'^F{uj) = C + 0{uj), and neglect linear 
and higher order powers of u. With this we obtain 

-ImE''P(w) =7rC|w| (17) 

which is linear in u for constant C. 

Similar formulas have been derived for the optical mass renormalization and scattering 
rate which when applied to the specific case of graphene give 

A^^c.) = ^ r dQa'F{Q)P r dJ-^^^^±^ In , + , (18) 
^ ' uj^Jo Jo 1^0 (u' + ny-iu^ ^ ' 



and in the limit a; — )■ we get back the same result as Eq. (14) namely X°^{ijj = 0) = A^^(a; 
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0). One should note that for finite Eq. (18) can be solved as a series of logarithms. The 
simplified approximate formulas derived above can be used to analyze the numerical results 
obtained from the evaluation of the full Kubo formula for the conductivity described in the 
previous section and to obtain insight into the nature of the effective excitation spectrum that 
scatters the electrons when electron-electron interactions are accounted for. The expression 
for the optical scattering rate is 



r°J'(a;) 



27r 
vr/io 







a^F{VL) 



LO — Vt 



1 + 



for Vt < u < -\- Vt 

for UJ > Hq + \ l. 



1^0 



(19) 



The zero u limit of 1/t°^{u) for the same model a^F{Q) as previously used gives 1/t°'^{uj) 
ttC\uj\ which is the same as its quasiparticle counterpart. 



IV. RECOVERING A BOSON SPECTRUM 



It is useful to use equation (19) in a maximum entropy inversion of our numerical data 
for the optical scattering rate which includes the electron-electron interaction shown in Fig. 
1 For a general kernel, K{u, Q), and input data, I{ou), with /(w) = K{u, Q)a'^F{fl)dfl 
the deconvolution of this equation to recover an effective spectral density, a'^F{Q) is ill 
conditioned and here we use a maximum entropy techniqud^. We begin by discretizing 
the equation to = K{i, j)a'^F{j)AQ where AQ is the differential increment on the 
integration over Qj = jAQ. We define a by 

1=1 * 

where lii) is the input data, S(i) = K{i, j)a'^F{j) the calculated value from the known 
kernel and a given choice of a'^F{Q), and at is the error assigned to the data Appro- 
priate constraints such as positive definiteness for the boson exchange function, a'^F{Q), are 
introduced and the entropy functional 

L = ^-aS (21) 
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FIG. 4. (Color online) The top frame displays the optical scattering rate 1/t°p{uj) (in meV) as a 
function of uj (in meV) for the case of pure EEI (solid black) . The bottom frame displays the various 
a'^F{Lo) obtained through maximum entropy inversion which are used to obtain fits to the input 
EEI of the top frame. The dashed blue curve is the maximum entropy fit to the solid curve while 
the dash-doted red curve (which does not fit to the data well) is obtained by removing spectral 
density at low uj, described in detail in the text, and shown in the lower frame as the dashed-dotted 
red curve assumed linear in uj below 100 meV. The solid green curve is the electron boson spectral 
density obtained from our maximum entropy fit. The inset in the top frame provides an expanded 
view of the scattering rates at small uj. 
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FIG. 5. (Color online) Upper frame, the quasiparticle scattering rate 1/t'^p{u) (in meV) vs u (in 
meV) for various cases. Negative and positive data refers to sign of frequency uj as the quasiparticle 
scattering rate can be different for positive and negative frequency cases. The dashed-dotted blue 
curve is the inversion fit to the negative energy EEI data; the dashed double-dotted green curve 
is the fit to the positive energy EEI data; and the red dashed curve is the inversion fit to the 
EEI+EPI case for negative energies. Lower frame, the electron boson spectral densities recovered 
from the maximum entropy fit for EEI+EPI (solid red), EEI negative frequencies (dashed dotted 
blue) and positive frequencies (dashed double-dotted green). 



where the Shannon- Jones entropy, is minimized with 



S 



a^F{VL) - m{n) - a^F{VL) In 



a''F{VL) 



(22) 
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FIG. 6. (Color online) Top frame, the quasiparticle scattering rate l/r'?P(a;) (in meV) vs oj (in meV) 
for EEI alone (solid black) , its maximum entropy fit (dashed blue) and the results of a modification 
in the electron boson spectral density (dash-dotted red). The inset provides an expanded view of 
the low uo region. Lower frame, the recovered electron boson spectral density solid blue curve 
through maximum entropy inversion of the data of the top frame. The dash-dotted red curve is 
the same but with an arbitrary modification i.e. linear dependence in uo from to 50 meV. 



The parameter a in equation (21) controls how close a fit to the data is obtained. The 



parameter m(f2) is here taken to be some constant value on the assumption that there 
is no a priori knowledge of the functional form of the electron-phonon spectral density 
a'^F{VL). While there is no guarantee that a boson exchange model can successfully reproduce 
consistently, quantitatively, and accurately all the details of our calculated data for the 
optical self energy in the case of the EEI application of the inversion process can still provide 
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useful insight into the excitation spectrum responsible for the scattering in this case. 



In Fig. |4]we show results obtained from the optical scattering rate l/r°^(ci;) of Eq. (10) 
which is based on the Kubo formula, Eq. ([s]), and the electron spectral density, Eq. (|5]), 
with quasiparticle self energy given by equation ([T]) which is based on a random phase 
approximation for the dynamical screening in graphene. A Kramers-Kronig transform is 
also implied. This is how the imaginary part of the conductivity is obtained from its real 
part given in equation (|3]). Both real and imaginary parts of the optical conductivity are 
required in constructing the optical scattering rate. In the top frame the input data on 
1/t°p{uj) in meV as a function of u up to 1000 meV is shown (heavy solid black). The 
maximum entropy fit for cr = 0.9 is shown (dashed blue) and is seen to be an excellent fit 
over the entire range of u considered. Here the chemical potential was set at fiQ = 1000 
meV for definiteness. The recovered spectral density a'^F{uj) is shown in the lower frame 
(solid blue). An important feature to notice is that below 100 meV a'^F{u) remains finite 
to rather low energies before going to zero at very low u (not seen on the scale used for this 
figure). This is characteristic of the Coulomb interaction and reminds one of the marginal 
Fermi liquid model for which a'^F{u) is constant for all u greater than a low energy cutoff on 
the order of the temperature. It is interesting to look, at the low u behavior more closely. 
The dashed dotted red curve shows a modified spectral density which has been made to 
vanish linearly below 100 meV. For this spectrum we obtain the dashed-dotted red curve 
for the scattering rate shown in the top frame. We see that this arbitrary modification 
has significantly altered our low energy fit and that to get agreement with the input data 
a reasonably constant value of the spectral density is indeed needed i.e. a marginal Fermi 
liquid like behavior is indicated. A reasonable first measure of the difference at small u in the 
two spectra of Fig. ^ (lower frame) is the value of the spectral lambda, A = 2 " ^^'^^ du, 
which emphasizes this region. In the best fit case A = 0.193 while in the other it is reduced 
to 0.135 but clearly the first of these values is preferred as the low u value of the data for 
1/t°^{uj) shows linear in u dependence in the range shown and this is taken as characteristic 
of coulomb interactions, as emphasized in the inset which shows a magnified view of the low 
u dependence only. 

A scattering rate taken to be linear in u at small u implies a constant value of a'^F{u) at 
small u. This is true for both optical and quasiparticle quantities. In Fig. ISlwe show results 



for the imaginary part of the quasiparticle self energy. As noted in equations (15) and (16) 
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the form of the integral determining this quantity in graphene depends on the sign of u. 
For finite chemical potential /iq, even the density of electronic states is not symmetric about 
a; = 0. In the top frame of Fig. |5]we show results for l/r'^^iuj) for both positive and negative 
uj data for EEI alone and for a combined case of EEI+EPI for the negative range only (heavy 
solid black curves). In all cases the maximum entropy fits (shown in the top frame) are very 
good. These are calculated from the recovered electron-boson spectral densities presented 
in the lower frame. We note first that all curves remain finite at small u on the scale of the 
figure. Of course the maximum entropy inversion respects the constraint that a'^F{ijj) = at 
a; = 0. As is increased the two pure EEI spectra show a broad maximum before dropping 
towards zero above uj ~ 1400 meV for the negative energy and around 800 meV for the 
positive case. The chemical potential is 1000 meV. The third spectrum, shown as the solid 
red curve, has the same EEI background as the pure electron case but has an additional 
sharp peak due to coupling to an Einstein phonon aX u = VLe = 200 meV. This EPI part 
is clearly seen over the EEI background. It is this EPI part which is responsible for the 
sharp rise in the 1/t'^^{uj) curve of the top frame for this case. Should the coupling be to a 
distribution of phonons rather to an Einstein mode then the rise aX uj = VLe would be more 
gradual, reflecting the details of the distributed phonon spectrum involved. It is clear, on 
comparison of Fig. |6] for quasiparticle self energies with the data in Fig. |4] for the optical 
case, that a boson exchange theory does not provide a perfectly consistent picture between 
these two data sets. Nevertheless the shape of the recovered spectra for EEI are quite similar 
and neither shows any evidence for important coupling to plasmon structure and both are 
more characteristic of coupling to a particle-hole continuum. In this regard it is important 
to note that the quasiparticle data inverted in Fig. |6] is for positive energy data at k — kp. 
In the upper frame of Fig. [6] we repeat our maximum entropy fit to the scattering rate data 
for positive energy but also add a second (red) dashed-dotted curve for a case where the 
a^F{uj) has been altered at small u to go linearly to zero below 50 meV rather than stay 
constant for the maximum entropy inverted curve (solid blue). This modification of the 
spectrum has changed the spectral A from 0.327 to 0.254. We note that these two values of 
mass enhancement recovered from quasiparticle data, while of the same order of magnitude 
as these recovered from optics, are definitely larger than their optical counterparts. Thus a 
single electron-boson spectrum does not fit simultaneously both quantities. It is important 
to recognize, however, that the shape of the functions involved are not different. We are 
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mainly noting a magnitude difference. In the inset of the top frame of Fig. |6] we show 
the deterioration in our fit to l/r'^P(a;) data that our arbitrary modification of the low uj 
behavior of the spectral density has caused, i.e. when the red dash-dotted curve for the 
modified a'^F{uj) is made linear below 50 meV with no other changes. We emphasize that 
once again the quasiparticle scattering rate shows a linear in u dependence at small energies, 
reminiscent of the marginal fermi liquid model characteristic of coulomb interactions in this 
case. 

In Fig. [7] we compare our results for the electron-boson spectral density obtained from 
positive oj quasiparticle self energy (solid red) with our optical conductivity data (dashed 
blue), for two values of chemical potentials. The top frame has /io = 400 meV and the 
middle flame is for /xq = 1000 meV. Comparing these two frames we conclude that in both 
the quasiparticle and optical cases the boson spectral density does not change its magnitude 
when the value of /xq is changed and that the horizontal energy scale u is itself proportional to 
/io- Also, as we have previously noted, the shapes of the extracted a^F{(jj) are not different, 
rather there is only a scale difference between quasiparticle and optics. In the lower frame 
of Fig. [7] we make this point more clearly. There the dash-dotted blue and solid red 
are respectively the optical spectral density and the averaged positive and negative energy 
quasiparticle spectral density with this latter quantity scaled to result in the same spectral 
A value. The curves are not appreciably different and indicate that in both quantities we 
are coupled to the same excitations. 

For Coulomb correlations the quasiparticle self energy depends on magnitude of the mo- 
mentum. In Fig. [8] we show results for the optical scattering rate when this momentum 
dependence is pinned to its value aX k = kp for all fc's in Eq. ([t]). As is clear from the 
top frame this approximation strongly influences the absolute magnitude of the scattering. 
Comparing the pinned (EEI, k = kp) and regular (EEI) cases we see that the curves do not 
differ much when one considers their dependence on u. Both can be fit to a boson exchange 
theory using maximum entropy methods (dashed double-dotted blue curve for EEI k = kp 
and dashed dotted red for EEI). Results for the recovered spectral density a^F{uj) are shown 
in the lower frame. The upper blue curve with \EEi,kp = 0.30 has very much the same shape 
as the lower red curve for EEI with XpEi = 0.144, although their A's differ by a factor of 2. 
We interpret this as strong evidence that in optics the boson spectrum of excitations that is 
sampled is that due to the electron a.t k = kp. By implication the plasmaron structures cor- 
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responding to k near zero, are not important as we have already conclude from consideration 
of the frequency dependence of the optical effective mass. Extracting the entire spectrum 
of excitations through inversion has allowed us to address this issue in more detail than has 
previously been possible. 

Another important point is made in Fig. |9} In the top frame we show results for the 
optical scattering rate 1/t°^{u) in meV as a function of u to 500 meV. The low solid black 
line is for EEI alone which the upper black line includes an electron phonon component 
(EEI+EPI). Note the sharp rise in 1/t°p{u) at a; = 200 meV which is the energy at which 
we have placed the Einstein phonon. The dashed red curve is our maximum entropy fit 
to the lower black curve and the dash-dotted blue to the upper one. In both curves the 
fits are very good over the entire energy range considered. The recovered electron-boson 
spectral densities are shown in the middle frame. The dashed red curve is for EEI alone 
and the solid blue includes both EEI and EPI. Note that the inversion procedure provides 
almost exactly the input value used for the mass renormalization coming from the phonons 
in that Xeei+epi — ^eei = 0.128 to be compared with the input value Xepi = 0.125. 
While the Einstein contribution is characteristically peaked about u = Qe = 200 meV, the 
EEI provides a corresponding smooth and low amplitude background extending to a; > 300 
meV on which the phonon peak sits. In the lower frame we show additional results for 
the recovered electron-boson spectral density in the same case but obtained from a second 
derivative techniqu^^^'^ rather than by maximum entropy inversiorP^. It is well known that 
in conventional metals in which electrons are coupled by the electron phonon interaction 
the second derivative of the scattering rate times u namely (^^^^^ j = W{u) is closely 
related to the input a'^F{u) in the energy range where it is non zero. Above the cutoff 
in a'^F{u)), W{uj) shows additional negative tails not part of the original spectrum. These 
results however depend on the assumption that the electronic density of states does not vary 
significantly on the energy scale involved in a'^F{uj). Nevertheless as seen in the lower frame 
of the figure this also holds reasonably well for graphene. Below roughly 320 meV, solid 
blue and dashed red curves for W{u) are qualitatively very similar to the solid blue and 
dashed red curves for shown in the middle frame. Both techniques reproduce well 

the main features of the underlying spectrum. Also, separate signature of both phonons and 
EEI are clearly seen although strictly speaking the spectra are superimposed. But, because 
they have quite distinct energy variations this allows for a reasonable separation of the two. 
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There are quantitative differences such as the cutoff around 300 meV occurring at lower 
energy in the maximum entropy inversion than in the second derivative technique. 



V. BOSON STRUCTURE IN INTERBAND CONDUCTIVITY 

So far we have discuss the intraband conductivity of graphene which is the part most 
closely related to the conductivity in ordinary one band metals. It is this piece which is most 
naturally analyzed in terms of an optical self energy and which has been found most useful 
to obtain information on the spectrum of excitation involved in the inelastic scattering of 
the charge carriers. But the conductivity in graphene also has a second component which 
can be nicely separated from its Drude component by increasing the chemical potential /ig 
(charging effects). The interband part in the free band case starts at photon energy Q = 2^q 
where it rises sharply to its universal background value cxo = 7re^/(2/i). When interactions 
are included the chemical potential shifts to higher energies for the electron-electron case 
and to lower energies for the electron-phonon and the onset broadens. Here we are interested 
in the EEL In this case the rising edge is rather broad as shown by the solid black curve of 



Fig. 10 labeled as interband conductivity and we note modulating structures around lo/ jj,Q ^ 
0.9 and 1.3 coming from the electron-electron correlations. These can be brought out more 
prominently by taking a first derivative of the interband conductivity which displays two 
peaks; a narrow, low energy peak and a broad, high energy peak. These two peaks correlate 
perfectly with peaks at negative energy in the electron spectral density for the lower Dirac 
cone at k = kp, A~{kF, —u). Transitions from the lower to upper cone at k = kp appear 
to dominate the onset of the interband transitions, but the electron-electron interactions 
clearly provide a substantial broadening of this absorption edge. 



VI. CONCLUSIONS 

We have calculated the effect of electron-electron (EEI) and electron-phonon (EPI) inter- 
actions on the intraband optical effective mass renormalization X°p{uj) and associated optical 
scattering rate 1/t°p{u) of graphene. These quantities are not independent but rather are 
related by a Kramers-Kronig transformation. They are encoded with information on the 
effective spectrum of excitations which scatter the massless Dirac fermions responsible for 
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charge transport. For the EEI the quasiparticle self energy is computed in hnear approxi- 
mation in the bare potential reduced by the substrate dielectric function (average over top 
and bottom dielectrics), which is then screened dynamically in random phase approxima- 
tion. For the EPI the corresponding self energy follows directly from a knowledge of the 
electron-phonon spectral density a'^F{Q) known in graphene from density functional theory. 
Once the electronic self energy is known the optical conductivity and optical self energy 
follows from a Kubo formula. For the simplest case of coupling to a single Einstein bo- 
son (Jlfi) mode it is found that the optical mass renormalization exhibits a characteristic 
easily recognize peak at a; = Qg. This simple observation serves as a qualitative guide for 
the determination of the distribution of phonons or bosons that may be involved in more 
complex cases. Examination of our numerical results for A°^(w) vs cu when only EEI are 
considered leads directly to the conclusion that the effective boson spectrum involved in 
the electron renormalization consists of a relatively uniform distribution extending over an 
energy range of the order of the chemical potential /i and scaling linearly with /i. Changing 
the value of // by application of a gate voltage in a field effect device or by some other 
means such as seeding donor or acceptor atoms on the surface of a graphene sheet, will alter 
linearly the range of the spectrum while at the same time leaving the value of the spectral 
A = 2 [a'^ F (u) / u!]du! unchanged. 

If in addition to EEI we include also an EPI piece characterized by an Einstein peak 
for largest contrast, we recover the original EEI background with superimposed, Einstein 
peak at a; = Qe- While the two distributions cannot be displayed separately, they can 
still be individually recognized because of their distinct distribution in energy. While we 
have used for best effect an Einstein peak, a distributed phonon spectrum would still be 
expected to have sharp peaks and stand out from the relatively unstructured electron- 
electron background. In addition the shape of the EEI part can be manipulated in a specified 
way through charging while one would expect the EPI part to be much less affected. 

The interband optical transition provide unstructured background to the optical response 
of graphene distinct and additional to its intraband Drude like response. This new contri- 
bution starts at photon energy ou — 2/1 in the bare band case. Correlations smear this rising 
edge and we find that EEI introduce modulating structures to this edge which are found to 
correlate remarkably closely with equivalent structures seen in the charge carrier spectral 
density —uj) for k — kp in the lower occupied valence band in the case considered. 
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Another important conclusion of our analysis is that the EEI structures identified in the 
intraband optical self energy is closely tied to the equivalent structure in the quasiparticle 
self energy for k equal or close to kp only. The plasmaronic structure most characteristic of 
the spectral function A{k,uj) for k around zero (i.e. around the Dirac point) will display no 
features in the intraband optical response. 
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FIG. 7. (Color online) Top frame, the maximum entropy results for the electron boson spectral 
density, a^F{Lo), obtained from positive energy quasiparticle data (solid red) compared with that 
obtained from optical data (dashed blue) for chemical potential /io = 400 meV. The middle frame 
is the same but for the case /io = 1000 meV. The lower frame compares optical and quasiparticle 
spectra for /io = 1000 meV with value of mass enhancement A for the quasiparticle case adjusted 
to be the same as for the optics. Here it is the average of positive and negative energy quasiparticle 
inversions that is employed. 
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FIG. 8. (Color online) Top frame shows the optical scattering rate 1/t°p{oj) (in meV) vs uo (in 
meV) for EEI (solid black) with maximum entropy fit (dash-dotted red). The lower frame shows 
the recovered electron boson spectral density, a^F{uj), for pure EEI (solid red) and the pinned 
k = kp case (solid blue). The solid blue is data for EEI but with the electron-electron self energy 
pinned at A; = /ci? for all values of k and the fit is the dashed double dotted blue curve. 
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FIG. 9. (Color online) Top frame shows the optical scattering rate 1/t°p{uj) (in meV) vs lo (in meV) 
for the case of combined EEI+EPI (solid black) with the maximum entropy fit to the numerical 
data (dash-dotted blue). This is to be compared with the pure EEI case (solid black) with its fit 
(red dashed curve). Middle frame gives the electron boson spectral densities recovered from the 
maximum entropy inversions. The solid blue is EEI+EPI while the dashed red is for EEI alone. 
This last curve provides a background above which there is an additional phonon peak at 200 meV 
for the combined case. The bottom frame is the same as middle frame but now showing W{uj), the 
second derivative of 2^ t°p(u;) ^i^h respect to uj, which is used to get an estimate of the underlying 
spectral density in the combined EEI+EPI (solid blue) and EEI (dashed red) cases. 



29 




FIG. 10. (Color online) The carrier spectral density A~{kF,uj) (solid green) for the lower Dirac 
cone at momentum k = kp as a function of normalized energy normalized by the chemical potential, 
uj/Ho- The solid black curve gives the interband contribution to the graphene conductivity and its 
first derivative is given by the dashed blue curve. 
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